Bias-Controlled Selective Excitation of Vibrational Modes in 
Molecular Junctions: A Route Towards Mode-Selective 

Chemistry 

Roie Volkovich^, Rainer Hartle^, Michael Thoss^, and Uri Peskin^ 
^Schulich Faculty of Chemistry and the Lise Meitner 
Center for Computational Quantum Chemistry, 
Technion- Israel Institute of Technology, Haifa 32000, Israel 
^ Institut fur Theoretische Physik und 
Interdisziplinares Zentrum fur Molekulare Materialien, 
Friedrich- Alexander- Universitat Erlangen-Niirnberg, 
Staudtstr. 7/B2, D-91058 Erlangen, Germany 

Abstract 

We show that individual vibrational modes in single- molecule junctions with asymmetric 
molecule-lead coupling can be selectively excited by applying an external bias voltage. Thereby, a 
non-statistical distribution of vibrational energy can be generated, that is, a mode with a high fre- 
quency can be stronger excited than a mode with a lower frequency. This is of particular interest in 
the context of mode-selective chemistry, where one aims to break specific (not necessarily the weak- 
est) chemical bond in a molecule. Such Mode-Selective Vibrational Excitation is demonstrated for 
two generic model systems representing asymmetric molecular junctions and/or Scanning Tunnel- 
ing Microscopy experiments. To this end, we employ two complementary theoretical approaches, 
a nonequilibrium Green's function approach and a master equation approach. The comparison of 
both methods reveals good agreement in describing resonant electron transport through a single- 
molecule contact, but also highlights the role of non-resonant transport processes, in particular 
co-tunneling and off-resonant electron-hole pair creation processes. 



PACS numbers: 73.23.-b,85.65.+h,62.25.Fg 



I. INTRODUCTION 



Mode-selective chemistry, that is the control of specific conformational changes or chem- 
ical reactions by directing energy into specific vibrational modes of a molecule, has been a 
major goal and challenge of modern chemical physics p~H3]. A variety of different routes to 
achieve mode-selective chemistry have been considered. For example, ultrafast excitation 
of molecules by laser pulses provided a route for directing energy transiently into specific 
mode-excitation or bond-cleavage in the presence of Intramolecular Vibrational Energy Re- 
distribution (IVR) processes • Nevertheless, on long time scales the statistical distribu- 
tion of energy usually prevails and the control of specific conformational changes or chemical 
reactions by directing energy into specific vibrational modes still represents a challenge [7j. 

In recent years, much interest has been devoted to the study of molecular systems out 
of equilibrium, in particular single-molecule junctions [8H23]. These junctions consist of a 
single molecule that is clamped between two metal or semi-conductor electrodes. If these 
electrodes are set to different electrochemical potentials by applying an external bias voltage, 
electrons tunnel from one lead to the other trough the molecule. The geometrical structure 
of the molecular bridge, due to its small size and mass, is very sensitive to charge fluctuations 
induced by these tunneling processes. The tunneling electrons thus strongly interact with 
the vibrational degrees of freedom of the junction [TU IT9l42~T| 1231427] . The distribution of 
vibrational energy on the molecular bridge, which results from these interactions, is highly 
correlated with the applied bias voltage and often deviates from a Boltzmann distribution 
[261 [281438] . In a single- molecule junction, it is thus possible to control a "non-statistical" 
non-equilibrium distribution of vibrational energy by an external potential bias. This offers 
an alternative route to mode-selective excitation (or mode-selective chemistry). Moreover, 
in the steady-state transport regime of a molecular junction, this could be achieved even 
without the cumbersome preparation of a specific initial-state. 

Employing generic model systems, we have recently shown [31] that in a single-molecule 
junction the excitation of vibrational modes can be, indeed, selectively controlled by the ex- 
ternal bias voltage. Considering molecules with an asymmetric orbital structure and vibronic 
coupling to specific normal modes, it was demonstrated that, by adjusting the external bias 
voltage, such a system can be driven into different non-equilibrium states with different lev- 
els of excitation for specific nuclear modes. Particularly, the excitation of a high frequency 
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mode can be tuned much higher than that of a low frequency mode, thus overcoming the sta- 
tistical distribution of vibrational energy that favors the excitation of low frequency modes. 
In this article, we extend our previous studies [31], where the principle of Mode-Selective 
Vibrational Excitation (MSVE) in a single-molecule junction was demonstrated for the first 
time. To this end, we outline in detail how the excitation of a single vibrational mode can be 
controlled by an external bias voltage. This bias-controlled excitation is then generalized to 
more than one vibrational mode, reviewing the basic MSVE phenomenon and emphasizing 
the role of intra-molecular interactions. In particular, we demonstrate MSVE in the pres- 
ence of electronically mediated mode-mode coupling [32] , which results from coupling of the 
vibrational modes to the same electronic state. These interactions induce energy transfer 
between the vibrational modes and tend to distribute current-induced vibrational excita- 
tion between the different modes. Similarly, electronic correlations, e.g. due to Coulomb 
repulsion, may influence MSVE by reorganizing the electronic population between specific 
electronic states at the molecule, which, to some extent and in different ways, is related 
to the phenomenon of MSVE. Extending our previous work, we investigate complementary 
models for asymmetric molecular junctions exhibiting MSVE. These models are related, for 
example, to experiments on single-molecule junctions performed with a Scanning Tunneling 
Microscope (STM) [131 H51 El QB1 123 ESHU] • 

To describe this nonequilibrium transport problem, we employ two complementary ap- 
proaches. The first one is a Master Equation approach (ME) that is based on the time- 
evolution of the reduced density matrix of the molecular bridge [2E1 EH I3J2 [33, H2H2E]. 
The respective equation of motion is evaluated strictly to second-order in the molecule-lead 
coupling, for which all resonant transport processes are included. Within such a frame- 
work, no approximations with respect to the interactions on the molecular bridge need 
to be invoked, in particular not with respect to electron-electron interactions or electronic- 
vibrational coupling. However, higher-order processes [321 EHHSS] like co-tunneling processes 
or the broadening of molecular levels due to the molecule-lead coupling are missing in this 
description. Note that master equation approaches that take into account such higher-order 
effects have already been put forward j^2l H5l I531455] . In the present work, however, we 
selected a Nonequilibrium Green's Function approach (NEGF) [211 I32H3H 1381 [661475] to 
account for the higher-order effects. Especially for the description of multiple vibrational 
modes [23J [321 El] , the NEGF methodology is typically more efficient. We employ a nonequi- 
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librium Green's function approach, which was originally proposed by Galperin et al. [67] 
and recently extended to account for multiple vibrational and multiple electronic degrees 
of freedom of the molecular bridge J32H2U EH]- Other theoretical approaches used to de- 
scribe electron transport through a single-molecule junction are based e.g. on scattering 
theory [761482] . path integrals [83H85] . multiconfigurational wave-function methods [86H88] . 
flux-correlation approaches [89] or exact diagonalization [90J. 

The paper is organized as follows. The Hamiltonian that we use to describe electron trans- 



port through a single-molecule junction is outlined in Sec. II A The ME and NEGF method- 
ologies that we employ to calculate steady-state observables of a biased single-molecule junc- 
tion are briefly described in Sees. II B and |II C respectively. The basic physical mechanisms 
for vibrational heating and cooling in an asymmetric molecular junction, which lead to 
MSVE, are reviewed and analyzed for a single-level conductor in Sec. HI A[ Thereby, the 
role of electron-hole pair creation processes, which constitutes an important cooling mecha- 
nism in a molecular junction [3H EH E7J, is discussed in detail. Moreover, the comparison 
of results obtained from our NEGF and ME schemes enables us to distinguish resonant 



and non-resonant contributions. In Sec. Ill B we discuss MSVE in two model systems, 
representing generic asymmetric molecular junctions: A single-molecule junction with an 
intrinsically asymmetric molecular bridge (model A) |34j, and a junction, where the bridg- 
ing molecule is asymmetrically coupled to the leads (model B). The latter scenario is typical 
for STM experiments. We demonstrate that the magnitude and the polarity of the external 
bias voltage can be used to direct vibrational energy into specific vibrational modes, even 
in the presence of intra-molecular interactions that tend to suppress the effect. Thus, we 
predict that in the steady-state transport regime of a molecular junction, a bias-controlled 
"non-statistical" distribution of vibrational energy can be realized, where modes with higher 
frequency (stronger bonds) can be much higher excited than modes with a lower frequency 
(weaker bonds). 
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II. THEORETICAL METHODOLOGY 



A. Model Hamiltonian 

We consider electron transport through a single molecule that is covalently bound to two 
metal leads. To this end, we employ a generic model Hamiltonian 

H = H c \ + -f/vib, (1) 

describing the electronic, H e \, and the vibrational degrees of freedom, H^, of this transport 
problem. 

The electronic part of H can be represented by a discrete set of electronic states, located 
at the molecular bridge (M), and by a continuum of electronic states describing the electron 
reservoirs of the left (L) and the right (R) electrode. Tunneling of electrons from one lead 
to the other is described by the following model Hamiltonian (h = 1): 

H c \ = ^2 € ™ a ln a m + ^2 U m,na ] m a maia n (2) 
m£M m<neM 

+ ^2^kb[b k + ^2 (yK,m^K,kb[a m + h.c). 

The energies of the electronic states in the leads, which are addressed by creation and 
annihilation operators b\ and bk, are denoted by e^. The energy of the mth electronic 
state located at the molecular bridge is given by e m . These states are populated with 
creation operators o) m and depopulated with annihilation operators a m . The molecular 
bridge is bilinearly coupled to the electrodes, ~ b\a m , with coupling strengths VK, m ^,K,k- 
For simplicity, these coupling strengths are factorized into an electrode term, ^x,k, which 
determines the electrode's (K=L,R) spectral density, J^r(e) = J2k&K \^K,k\ 2 ^( e ~ e k)i an d 
a molecular term, which represents the coupling of the mth molecular state to the i^th 
electrode, vx, n - To model the leads we use a semi-elliptic conduction-band with a band- 
width of 47 such that the corresponding level-width functions read: 

£ 2 / 2 

rV, m „(e) = 1-nv K)m v K ^J K [e) = v K ^ m v K ^ — \J^ 2 - (e - fi K ) 2 . (3) 

Charging energies, e.g. due to Coulomb interactions, are accounted for by Hubbard-like 
electron-electron interaction terms, Um^a^amO^an. The Fermi-energy of the overall system 
is given by e F = eV. 
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Vibrational degrees of freedom of the molecule are described as harmonic oscillators, 

v m£M;i/ 

+ up v d\ v dp v + Y Vv,0„( d l„ + rf ftJ(4 + c ^)> ( 4 ) 

where the ladder operators c\,/c u address the ^th vibrational (normal) mode of the molecular 
bridge with frequency Q u . Changes in the nuclear potential energy surface due to electronic 
transitions between the single particle states are assumed to be linear in both the vibrational 
displacements, q u = ^(c^ + cj,), and the electronic densities, (^ m a m . The respective coupling 
strengths are given by \ v ^ m - To incorporate vibrational relaxation effects in a phenomeno- 
logical way [9T| l92]. each intramolecular vibrational mode is coupled to a thermal bath. The 
creation and annihilation operator for a bath mode with frequency up v are denoted by eft 
and dp v , respectively. The corresponding mode-bath coupling strengths are given by r\ v ^ v . 
All properties of the bath that influence the dynamics of the system are determined by the 
spectral densities J v {oj) = ^2 a vl^i 00 ~ 111 the calculations presented below, we use 
an Ohmic bath model with a cutoff frequency, u cv , and 



J v {u) = ^-u)e- u ' u °>». (5) 



B. Reduced Density Matrix Approach 

Weak coupling of the molecule to the electron reservoirs (the leads), as well as to the 
energy reservoirs (the nuclear bath modes), is essential in order to relate the junction's 
steady-state observables to intrinsic properties of the molecular bridge |93j . For this pur- 
pose it is instructive to regroup the different terms in the Hamiltonian into the molecular 
"system", which include the molecular electronic states and vibrational modes, and "bath" 
terms, which include the electron reservoirs and the nuclear baths. The full Hamiltonian is 
therefore rewritten as: 

H = H s + H B + # SB , (6) 
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with 

meM m<n£M ^ mgM;i/ 

#sb = X^,/3,(4, + + c *0 + X (v K , m ^K,k-bla m + h.c.) (9) 

XeL,R;fceX;meM 

The steady-state of the system under bias can be calculated by following the time-evolution 
of the system to its stationary state, starting from an arbitrary initial state. We consider 
an initial density matrix in a product form, 

p(0) = Ps(0)®/3 B (0), (10) 
Pb(0) = J] Pb,k®\[pb, v , (11) 

where ps(0) is any normalized system density matrix with trg[ps(0)] = 1, and trs[...] de- 
notes the trace over the subspace of the molecular conductor. The electronic reservoirs are 
described by a product of equilibrium density operators, 

P B ' K = r - ^ — ; — , , h t h ' ( 12 ) 

tr[e~ ^ fceK {ek ^ fiK > b k bk ^ 

with K G L,R. Thus, we assume the leads to be in a thermal equilibrium state, which is 
characterized by the temperature k^T and the electrochemical potentials fix- Similarly, the 
nuclear baths are represented by a product of equilibrium density operators, 

Pb,» = r^; — ,t . » ( 13 ) 

tr [e _ ^ ^- 3 " Pv ^ ^ ] 

where v denotes the molecular mode to which the particular bath is coupled. The exact 
time-evolution of the full density operator is given by the Liouville-von Neumann equation, 

|p(t) = p(t)]. (14) 

Transforming the respective operators to the interaction representation, 

6\t) = e i[Hs+H B ]tQ^ e -i[H s+ H B ]t^ (15) 
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the Liouville-von Neumann equation can be rearranged to give [56] 



dt 



p\t) = -iiHUt),?®] - f dt'lHUt), [HUt'lp 1 

Jo 



—i 



f dt' f df[Hl B (t), [Hl B (t% [HUt'% p\t")]]]. (16) 
Jo Jt' 



Assuming weak system-bath coupling, the third term in Eq. (16) can be neglected, which 
yields 

-p\t) = -i[Hl B (t),p\t)}- J dt'[Hl B (t),[Hl B (t'),p\t)]]. (17) 

Defining the reduced density operator p\(i) = ti B [p l (t)], one obtains the well-established 
(Markovian) Master equation for p s (t) [351 S3 El EH EM] by replacing p l {t) by p B (0) <g> 
Pg(t) in Eq. (17), and taking the integration limit to infinity, /J — > J °°, 

POO 

-p s (t) = -t[H s , Ps (t)]- j dt>ti B [H S B,[H SB (t-t'),p B (0)p s (t)]], (18) 



d_ 

Ft' 



with H sb (t) = e -i(Us+HB)T H SBe i(H s +H B )T _ Thereby, we used that [p B (0),H B ] = and 
^b[pb(0)Hsb] = [56J. To evaluate this equation of motion, it is convenient to use the 
eigenstates of the molecular system Hamiltonian, 

H s \l)=E l \l). (19) 



Taken in this basis, Eq. (18) corresponds to the Redfield equation JDUISSl 97J. The molecular 
system observables at steady-state can be calculated from the infinite time limit of ps(t)- 
Moreover, the effects of coherences between the system eigenstates can be neglected in this 
limit as long as the molecular levels are non-degenerate [35] (Coherences between quasi- 
degenerate molecular levels can play an important role, as in the case for molecular motors 



[57]). For the model systems that we study in Sec. Ill coherences are not important in the 
steady state limit and consequently, the equation of motion for the diagonal matrix elements 
of the reduced density matrix, that is the populations Piif) = ps,i,i(t), is given by 

JU(t) = $> R + n h + J2 * {v) h>Pv{t)- (20) 

V v 

The respective rate matrices for electron tunneling take the form, 

[^ L/R k, = (i - ^MT h + r K R;e ) - <v E rtf* + r M /R ;e )> (2i) 



and are determined by the spectral densities and the Fermi occupation numbers at each 
electrode, 

ltf c/h ee 2vr| ^t; L/R , m [ay^| 2 J L/R (^ - E l )f^{E v - E t ), (22) 

m 

with 

Ji /R (e) = 1-/ C L/R (e). (24) 

Similarly, the rate matrices describing the coupling of the vibrational modes to their thermal 
bath take the form, 

= (i - M(rO? + rj£) - 5 hl , (1$ + rjj*), (25) 

and are accordingly determined by the spectral densities and the phonon occupation numbers 
for each system mode, 

T^ /d = 2tt|[4 + c^tfj^Ev - Ei)n*/ A {E V - E t ), (26) 

with 

n d = (e) = n u (e) + 1. (28) 
Observables of interest, such as the steady-state current from left to right, 

/ L ^ R = lim ^^2e[^,]P^K ( 29 ) 

Z J' 

the average level of excitation of mode z/, 

(4c,) = iimVp^)(/:|4 Ci ,|0, (30) 

and the populations of the electronic states, 

(c4a m ) = lim y^Pi(t){l\a) m a m \l), (31) 



z 



are calculated from the infinite time limit of the Pi(t). Thereby, ni is given by n\ = 

J2 m ( l \ a tn a m\l)- 



C. Nonequilibrium Green's Function Approach 



Alternatively, vibrationally coupled electron transport through a single-molecule junction 
can be described employing a nonequilibrium Green's function approach. Using such an 
approach facilitates the description of higher-order effects by the associated Dyson-Keldysh 
equations. The comparison of results obtained from the reduced density matrix approach and 
the nonequilibrium Green's function approach allows to elucidate the role of these effects in 
vibrationally coupled electron transport through a single-molecule junction. Here, we apply 
the method originally proposed by Galperin et al. [HZ], which we have recently extended 
to account for multiple vibrational modes and multiple electronic states J32HHI1 EH]- The 
approach is based on the small polaron transformation of the Hamiltonian H [28j E21 EH] 

H = e s He' s (32) 

= ^ ^rn a l n a m 4" ^ ^u C u C u + ^ U m ^CL^Q, m^rfln (33) 
m v n<m 

+ Yl ekb k bk + Yl w ^4a 

k v,f3 v 

+ Y ( V K,rr£K,kX m b{a m + h.C.) + ^ ^M\ u + ^J(4 + C v) , 



fceL,R;m£M vfi v 



with 



S = -l Y ^ a ln a mP V , (34) 



mu 



X m = exp[i^2~^-p v ], (35) 

v v 

Vv = -4= {cu - 4) • (36) 

The transformed Hamiltonian, H, thus contains no direct electronic-vibrational coupling 
term, but polaron shifted state-energies e m = e m — X/i/C^mAV)' vibrationally induced 
electron-electron interactions, U m ,n = U m , n — 2 ^^A^A^/^), and shift operators X m 



that renormalize the molecule-lead coupling term. Note that in Eq. (32) we have neglected 
the renormalization of the molecule-lead coupling term due to coupling of the vibrational 
modes to the thermal baths (67]. Furthermore, the renormalization of the electron-electron 
interaction terms, Um^a^ama^an, due to these interactions are also discarded. Such bath- 
induced renormalizations are beyond the scope of this paper. Also note that the small po- 
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laron transformation does not allow for an arbitrarily strong coupling between the vibrational 
and the bath modes [67J, which means for the given spectral densities that Q, < VL v u c ^/A. 

The single-particle Green's function G m>m i(T, r') is the central quantity of (nonequilib- 
rium) Green's function theory. With this Green's function all single-particle observables, 
e.g. the population of levels or the current through a single- molecule junction, can be read- 
ily calculated. For the computation of the single-particle Green's function G mim i{r,r') we 
employ the following ansatz [32H3^ [38| 167] : 

G m>m '{r,r') = -i(T c a m (r)aJ n ,(r / ))^ (37) 
= ^(T c a m (T)XUr)al(T')xUr')) R (38) 
w G m , m/ (r,r')(T c X m (r)Xl l (T')) B , (39) 

with the electronic Green's function G mjm /(r, r') = — ?(T c a m (r)a] n ,(r / )) i 5- and T c the time- 
ordering operator on the Keldysh contour. The indices H/H indicate the Hamiltonian, 
which is used to evaluate the respective expectation values. The factorization of the Green's 
function G m ^ m i into a product of an electronic correlation function, G mjm >, and a correlation 
function of shift operators, (T c X m (r)X r ' rt/ (r / ))jj, is justified, if the dynamics of the electronic 
and the vibrational degrees of freedom are decoupled. This is conceptually similar to the 
Born-Oppenheimer approximation |100t 110 J] . Accordingly, for transport through a single- 
molecule junction, an (anti-)adiabatic regime is defined by T K)mm 3> Q (T K ^ mm <C fl). 

The self-energy matrices for the electronic part of the Green's function can be determined 
from the equation of motion 

(id T - e m )G mtm >(T, T')(-id T > - e m ) = S(t, T')(-id T > - e m ) (40) 

+ ^L,my(l", t') + ^R,m,m'(T, t') + Scoul,m,m' ( r ; T ')- 

Here, self-energy contributions due to the coupling of the molecule to the left and the right 
leads are denoted by £l )TO , to '(t, t') and SR, m>m '(T, r'), while correlations that result from the 
electron-electron interaction term, Um^a^ama^an, are summarized in £coui,m,m/ ( r 5 t 1 ). We 
treat this latter part of the self-energy in terms of the elastic co-tunneling approximation [33l 
EH1 [Z2 H02]. We therefore approximate S Couljm , m /(r, r') by the self-energy S^ oulm m , (r, r'), 
which describes electron-electron interactions in the isolated molecule exactly. The self- 



11 



energies that describe the coupling of the molecular bridge to the leads, 

S L /R,m,m'( r > r/ ) = E t; L/R,m^L/R,m'|^L/R,fc| 2 fl'fc(r,r / )(T c X m /(r / )X^(r)) i? , (41) 
fceL/R 

are evaluated up to second order in the molecule- lead coupling, where ^(r, r') denotes 
the free Green's function associated with lead state k. The real-time projections of these 
self-energies determine the electronic part of the single-particle Green's function. In the 
energy-domain the corresponding Dyson-Keldysh equations read 

Gtm^) = G°S(t) + E G °S(t) (^(6) + S£, n ,(6)) G^ m ,(e), (42) 

n,n' 

Gl%(e) = J2 Gr m,n(t) (sf&to + £r^)) ^, m ,(e), (43) 



with 



(44) 



G° m t(e) = g^M + E*n a ( e ) S clt^(^)^(6), (45) 



n,n' 



ti» - 4^— W - (46) 



For the computation of 

n m ) 1_(p « )m n^ )m ) ^ rr ^ ^ (47) 

we use the populations of the electronic levels 

Im [G< m (r = 0)] , (48) 



G^(e) = W E flK 1 

ct=1..2 dim ( M ) VmSM 



n,, 



that we determine self-consistently, and vectors p a that point to the edges of a dim(M)- 
dimensional unit cube. 



Eqs. (42) and (43) give the exact result in the non-interacting limit, where X UjTn — > and 
U m ,n 0. One should bear in mind, however, that the elastic co-tunneling approximation 
treats the eigenstates of H effectively as independent transport channels. This description 
therefore needs to be applied with care, if these channels cannot be treated independently 



from each other, e.g. in the presence of quantum interference effects [38 1I103H108] . Moreover, 
Kondo physics [59l462| 165] is also not included in this description, as it employs a (self- 
consistent) second-order expansion in the molecule-lead couplings. 

12 



The correlation function of the shift operators is obtained using a second-order cumulant 
expansion in the dimensionless coupling parameters j^™ [321 EH EZ] 

(WAM)i = -p (^^f^y^ ~ i K ^nJt'* DvAt,t) ) ' (49) 

V v,u' V V V V / 

where we use the momentum correlation functions 

D u y = -i(T c p u (T)p u/ (T')) s . (50) 
Employing the equation of motion for D vv i 

' d 2 T - nl)D v y(T,T')(-d 2 T , - Hi,) = 6(T,T>)(-d 2 T , - n» (51) 



+n ba tvy (t, r') + Ilei^^^r, t'), 

we determine the corresponding self-energy matrices n bathjVjl/ / and U e i,uy- The self-energy 
matrix 

nbath^y = ^^kjd^lM (52) 

ft 

includes the coupling of the vibrational modes to the thermal baths, where D°o denotes the 
free Green's function of bath-mode (5 U . The electronic self-energy part, U e [ iU y, describing 
the interactions between the vibrational modes and the electronic degrees of freedom of the 
molecular bridge, is evaluated to second order in the molecule-lead coupling [321 EH E7J 

n e i )IV / (t,t') = -i ^ ^ u ^" ,m (E m;m /(r, r')G m/>m (r',T) + E m ' !Tn (r', r)G mim '(r, t')). (53) 



m,m' 



Thereby, we use the noncrossing approximation, where contributions mixing mode-bath 
and molecule-lead couplings are disregarded. Since Ii e \,u,u' depends on the electronic self- 
energies S mjm / = SL, m ,m' + £R,m,m' and Green's functions G mtm ', the respective Dyson- 
Keldysh equations need to be solved iteratively in a self-consistent scheme [32| [67] . 

With these Green's function, D u y and G m ,m', the vibrational excitation of each vibra- 
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tional mode is obtained according to 



(a v + Im [D< u (t = 0)] - (b v + (54) 



+E^ Im [G<, m (t = o)] 



m<m' 



i>. 



+2 £ ^MG< m (t = 0)]Im[^(t = 0)] 



"2 £ ^M^>(^0)]Im[G<, m ^O)], 

m<m' 

= E^riN^( 1 + 2 ^(^))' ( 56 ) 

which is consistent with the second order expansion used for the evaluation of the associated 
Green's functions. In Eq. (54), we also use the Hartree-Fock factorization: {a i m a m a [ m ,ami) ~ 



{(^ m a r ^{a' m ia m i) — (aj n a m ')(a] n ,a m ). The current is calculated employing the Meir-Wingreen 
formula [109J 

1 = 2e / ^ E ( S L, m , m '(<0G>,, m (e) - S> mim ,(e)G<, m (e)) . (57) 

m,m' 

Note that this scheme, including the elastic co-tunneling approximation, is current- 
conserving. 



III. RESULTS 



In this section, we employ the ME and the NEGF methodology, outlined in Sees. II B 



and |II C[ respectively, to analyze transport characteristics of asymmetric single-molecule 
junctions. In particular, we consider a single vibrational mode coupled to a single electronic 
state in Sec. |III A and demonstrate that an external bias voltage can be used to control its 



level of excitation. In Sec. Ill B , we extend these studies to two vibrational modes and show 
that their level of excitation, despite their different frequencies, can be selectively controlled 
by switching the polarity of the applied bias voltage. To demonstrate the generality of the 
phenomenon and to corroborate our findings, we consider a variety of parameter regimes, 
including the effect of intra-molecular interactions. 
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TABLE I: Model Parameters (Energy values are given in eV, K £ {L,R}) 
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A. Bias-Controlled Excitation of a Single Vibrational Mode in an Asymmetric 
Molecular Junction 

We first consider a model with a single vibrational mode coupled to a single electronic 



state in Sec. Ill A and study how the excitation of the vibrational mode can be controlled 
by an external bias voltage $. To this end, we consider a vibrational mode with frequency 
Q = 0.15 eV that is coupled to a single electronic state with coupling strength Ai t i = 0.6Q\. 
The electronic state is located e\ = 0.6 eV above the Fermi-level, and is asymmetrically 
coupled to the left, -Ul,i = 0.1, and to the right lead, v^i = 0.03, respectively. The specific 
model parameters are detailed in Table [TJ These parameters (including those given in Tab. 
Til) reflect typical values for molecular junctions, as they are determined for example in 



ab-initio calculations [801 fTT0HTT2| ITT2HTT5] or experiments [H CEHrEEJ E3H27]. 



Current- voltage characteristics for this model molecular junction and the respective pop- 
ulation of the electronic state are shown in Fig. [TJ Thereby, the solid black lines represent 
results obtained with the reduced density matrix approach, while the dashed black lines 
show results, for which we employed the nonequilibrium Green's function approach. The 
results of both approaches agree very well. Minor deviations between the approaches occur 
due to the broadening of the molecular levels, which is not included in the ME scheme. For 
positive bias voltages, the current and the population of the electronic state display a single 
step. This step indicates the onset of transport at e$ = 2ei (ei = t\ — Af denotes 



the polaron-shifted energy of state 1 (cf. Sec. II C)), where electrons from the left lead can 



resonantly tunnel onto the molecular bridge. Notice that the molecular energy level is lo- 
cated well above the Fermi-level of the junction, that is by several units of the vibrational 
frequency: Z\ — ep > 30i. Therefore, at the onset of transport by electron tunneling from 
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the left electrode, several inelastic channels corresponding to processes described in Figs. 
[2^i, and [2]i open up simultaneously. As the bias increases further, additional heating 
channels become available, involving tunneling of high energy electrons from the strongly 
coupled (left) electrode onto the molecule (see Fig. [2Jd) . However, these additional channels 
do not significantly increase the current, since the bottleneck for transport in this asymmet- 
ric junction is tunneling processes from the molecular bridge to the weakly coupled (right) 
electrode that are already active. Accordingly, the electronic state is populated (from the 
left) much faster than it is depopulated by tunneling processes to the right, and is therefore 
almost fully occupied for e$ > 1e\. For negative bias voltages, however, both the current 
and the electronic population exhibit a number of pronounced steps at $ = — 2(ei + nfl\) 
(n G No). Again, different inelastic transport channels open up simultaneously at the onset 
of the current. However, in this the bias decreases further, $ < — 2ei, additional 

heating channels open up one by one at the bottleneck for transport, that is additional tun- 
neling processes with respect to the right lead. Since these processes are inactive for higher 
negative bias voltages, — 2ei < $ < 0, one observes significant steps in the current-voltage 
as well as in the respective population characteristics. Notice that in this bias direction the 
electronic state is depopulated (to the left) much faster than it is populated (from the right), 
so that it is almost unoccupied for e$ < — 2ej.. The relative step heights that occur in these 
characteristics qualitatively reflect the transition probabilities ^ ( -$^M e - ^ 1,1 / 01 ' for a 
transition from the vibrational ground- to its nth excited state. For a quantitative analysis 
of the step heights, however, the nonequilibrium state of the vibrational mode, which is typ- 
ically highly excited, needs to be considered (cf. Fig. pi). As a result of electronic- vibrational 
coupling and the asymmetry in the coupling to the leads, the current-voltage characteris- 
tics thus exhibits a pronounced asymmetry with respect to the polarity of the applied bias 
voltage $, which is also referred to as vibrational rectification. This has been theoretically 
analyzed [33, EE] and experimentally verified [HI [T71 [T9] before. 

In contrast to these electronic observables, the corresponding vibrational excitation num- 
ber, (cl,c v )jj, increases in a series of distinct steps for both polarities of the bias voltage 
(see Fig. [3]). This finding cannot be solely understood in terms of electron transport pro- 
cesses. Although vibrational excitation is a result of inelastic electron transport processes 
(cf. Figs. [2k-d), another class of processes, which does not contribute to the current, needs 
to be considered. Resonant electron-hole pair creation processes [3U ESI EZ] , such as those 
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FIG. 1: Upper Panel: Current-voltage characteristics for a model molecular junction comprising a 
single electronic state that is coupled to a single vibrational mode and asymmetrically to a left and 
a right lead. The solid black line is obtained with the ME approach, while for the dashed black line 
NEGF is used. Lower Panel: The corresponding population of the electronic state as a function 
of the applied bias voltage <£. The asymmetry of the population characteristics with respect to 
the polarity of the bias voltage is a result of the asymmetric molecule-lead coupling. Electronic- 
vibrational coupling translates this asymmetry also to the respective current- voltage characteristics, 
which otherwise (i.e. without electronic- vibrational coupling) would be almost anti-symmetric with 
respect to 
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FIG. 2: (Color online) Basic schemes of vibrationally coupled electron transport through a single 
electronic state. Panels a) and b) depict examples for emission processes, where an electron se- 
quentially tunnels from the left lead onto the molecule and further to the right lead, thereby singly 
exciting the vibrational mode of the molecular bridge (red wiggly lines). Such emission processes 
are effectively 'heating' the junction (local heating). Examples for respective absorption processes 
are shown in Panel c) and d), where electrons tunnel from the left to the right lead by absorbing 
a quantum of vibrational energy (blue wiggly line). 

depicted by Figs. |4^i and [4}d, are effectively cooling the vibrational mode and diminish the 
current-induced vibrational excitation. Since these processes involve two sequential tunnel- 
ing events, they occur with the same probability as respective transport processes. Due to 
the asymmetry in the molecule-lead coupling, electron-hole pair creation processes with re- 
spect to the left lead are the most important ones. They are typically more effective the less 
vibrational quanta are involved. For this particular model system, cooling by electron-hole 
pair creation is therefore more pronounced for positive bias voltages, $ > 2e 1; where e.g. an 
electron-hole pair in the left lead can be produced by absorption of just a single quantum 
of vibrational energy. For negative bias voltages, $ < — 2e 1; the creation of an electron- hole 
pair in the left lead requires the absorption of more than ten vibrational quanta. As a result, 
vibrational excitation is much smaller for positive bias voltages than for negative ones. In- 
creasing the bias voltage, these pair creation processes are blocked one by one, as the energy 
gap between the molecular level and the electrode chemical potential increases. The steps 
in the vibrational excitation characteristics thus become larger with increasing bias voltage 
due to less efficient cooling by electron-hole pair creation processes [351 3TJ . This blocking 
of pair creation processes appears for both polarities of the bias voltage. 
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FIG. 3: (Color online) Vibrational excitation characteristics corresponding to the current- volt age 
and population characteristics shown in Fig. [TJ Due to the asymmetry in the molecule-lead coupling 
of this model molecular junction, electron-hole pair creation processes are cooling the vibrational 
mode more efficiently for positive bias voltages than for negative ones, leading to a strongly asym- 
metric excitation characteristics. The external bias voltage thus can be used to control the level of 
vibrational excitation in this asymmetric molecular junction. 

Apart from the broadening of steps, the ME and the NEGF approach give almost the same 
vibrational excitation characteristics. However, NEGF can be expected to give slightly larger 
values for the vibrational excitation, because inelastic co-tunneling processes [32 | l33~ t |53 | 163] . 
which are not included in the ME scheme, additionally contribute to the level of excitation for 
the vibrational mode. This is particularly important in the off-resonant transport regime, i. e. 
for | $ | < €i, where NEGF gives a small vibrational excitation while ME does not. Significant 
deviations between both approaches occur only for large bias voltages. Especially for large 
positive bias voltages, e.g. at $ > 2V, the vibrational excitation obtained from NEGF is 
significantly smaller than the one obtained by the ME scheme. We attribute this behavior 
to the contribution of cooling by off-resonant electron-hole pair creation processes, which 
are missing in the ME approach (an example of an off-resonant pair creation process is 
depicted in Fig. Ep). These processes become the dominant cooling mechanism at large bias 
voltages, where resonant electron-hole pair creation processes are suppressed, as they require 
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FIG. 4: (Color online) Example electron-hole pair creation processes in a molecular junction. Panel 
a) depicts an electron-hole pair creation process with respect to the left lead by absorption of a 
single vibrational quantum. Panel b) represents an electron-hole pair creation process with respect 
to the right lead by absorption of two vibrational quanta. The absorption of two vibrational quanta 
typically occurs with lower probability. Panel c) shows an off-resonant electron-hole pair creation 
process. 

increasingly higher vibrational energy. In contrast, off-resonant pair creation processes can 
occur by the absorption of just a single quantum of vibrational energy for all bias voltages. 
The ME approach thus gives a somewhat larger vibrational excitation than the NEGF 
method for bias voltages, where resonant electron-hole pair creation processes are strongly 
suppressed. 

The importance of cooling by electron-hole pair creation processes can be corroborated by 
studying the behavior of the vibrational excitation characteristics with respect to the energy 
of the electronic state, which for a given bias voltage influences the efficiency of resonant 
electron-hole pair creation processes. Fig. [5^, shows results for the vibrational excitation 
as a function of the energy e% at a fixed bias voltage. Thereby, the blue lines refer to a 
fixed bias voltage of $ = — 2V, and the red lines to $ = 2V. As before, solid (dashed) 
lines refer to calculations performed with the ME (NEGF) scheme. If the energy of the 
electronic state is closer to the Fermi-level of the system, ei — > 0, resonant electron- hole 
pair creation processes are more strongly suppressed, as they require the absorption of an 
increasing number of vibrational quanta. The less efficient cooling by resonant electron-hole 
pair creation leads to the general observed trend of an increasing vibrational excitation with 
a decreasing energy of the electronic state (for example from ei = leV to t\ = 0.4 eV). For 
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t\ < 0.4 eV, the NEGF scheme gives a significantly smaller vibrational excitation than the 
ME method. Here, cooling by off-resonant electron-hole pair creation processes (missing 
in the ME treatment) results in a significantly lower level of vibrational excitation. Notice 
that for an asymmetric junction, the value of Z\ for which the vibrational excitation obtains 
its maximal value differs from zero and depends on the bias polarity. Considering, e.g. 
the negative bias voltage (blue lines), a maximum of vibrational excitation is obtained 
at e"i ~ 0.3 eV (ei ~ 0.35 eV). For this bias, shifting the electronic level to lower values 
enhances cooling at the left electrode and suppresses cooling at the right electrode. Due 
to the asymmetry in the molecule-lead coupling, pair creation processes with respect to the 
left lead are more important than with respect to the right lead. A minimum of cooling 
efficiency by electron-hole pair creation processes, which corresponds to a maximum in 
vibrational excitation, is thus reached for positive values of Z\. Similarly, for positive bias, a 
maximum in vibrational excitation appears for negative values of Z\. Fig. [5]d represents the 
ratio (c{ci)$ = _2v/( c i c i)<i>=+2V- It shows that the asymmetry in the vibrational excitation 
characteristics, as well as in the current- volt age characteristics and the electronic population 
(data not shown), disappears, once the electronic level is located close to the Fermi-level of 
the system. This demonstrates that for an electronic state close to the Fermi-level, which can 
be controlled for example by a gate electrode [201 IH6H118] , the efficiency of both current- 
induced heating and cooling by electron-hole pair creation processes is the same for both 
polarities of the bias voltage $. 

At this point, it is interesting to study the extent of vibrational excitation for different 
ratios of the molecule-lead couplings t^i/t^i. In Fig. [6j we show the level of excitation 
of the vibrational mode as a function of the ratio vr,i/vl,i, where = 0.1 eV is fixed. 
Again, red and blue lines refer to calculations performed at a fixed bias voltage of $ = ±2 V, 
respectively. Trivially, for a symmetric junction with v^i = Vl,i> we obtain the same level of 
vibrational excitation for both polarities of the bias voltage. Decreasing the coupling to the 
right lead, the asymmetry in vibrational excitation increases almost linearly. Interestingly, 
for negative bias voltages, the vibrational excitation obtained by the NEGF scheme (dashed 
blue line) and by the ME method (solid blue line) approach one another upon decreasing 
■ur 5 i. This points to the fact that off-resonant electron- hole pair creation processes with 
respect to the right lead become strongly suppressed. For even smaller coupling strengths 
to the right lead, v^i < 0.03 eV, the turnover in the dashed blue line (NEGF scheme) 
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FIG. 5: (Color online) Upper Panel: Vibrational excitation as a function of the energy, e\, for 
a molecular junction with an electronic state coupled to a single vibrational mode and asym- 
metrically to the leads. Red and blue lines refer to results that are obtained for a fixed bias 
voltage, $ = ±2 V, respectively. Lower Panel: Ratio of the average vibrational excitation numbers 
(c|ci)$=_2 v/( c i c i)*=+2 v shown in the upper panel. The further the electronic level is located 
from the Fermi-level of the system, which can be controlled for example by a gate electrode, the 
more pronounced is the asymmetry in vibrational excitation for the different polarities of the bias 
voltage $. 
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FIG. 6: (Color online) Vibrational excitation of a molecular junction with a single electronic state 
coupled to a single vibrational mode for different coupling strengths to the right lead, vn t \. Thereby, 
the coupling to the left lead is fixed, = 0.1 eV, as is the bias voltage, <3? = ±2 V, for the red and 
the blue lines, respectively. Asymmetric coupling to the leads facilitates the control of vibrational 
excitation by the polarity of an external bias voltage. 

indicates that in the limit v^i — >■ current-induced vibrational excitation vanishes, as does 
the corresponding current. The solid blue line (ME scheme) exhibits the same turn-over but 
for even smaller coupling strengths to the right lead. On the other hand, for $ = 2V, the 
dashed and solid red lines remain well separated upon decreasing the coupling to the right 
lead, as the ratio between resonant and off-resonant pair creation processes with respect to 
the left lead remains constant. 

So far, we have discussed cooling mechanisms for the vibrational mode, which are solely 
induced by electronic-vibrational coupling. Other vibrational energy relaxation processes, 
which can be of relevance in molecular junctions, include Intramolecular Vibrational En- 
ergy Redistribution processes (IVR) or energy transfer to the environment (e.g. phononic 
excitation of the electrodes) [131 |82| I119r fl21j. Such relaxation mechanisms are commonly 
described by coupling of the primary vibrational mode(s) to a thermal heat bath. Fig. [7] 
shows the level of vibrational excitation as a function of the mode-bath coupling strength 
Ci for a fixed bias voltage (blue lines correspond to $ = — 2 V, and red lines to $ = +2V). 
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FIG. 7: (Color online) Vibrational excitation as a function of the mode-bath coupling strength, 
(i, for a molecular junction with an electronic state coupled to a single vibrational mode and 
asymmetrically to the leads. Red and blue lines refer to results that are obtained for a fixed bias 
voltage, = ±2V, respectively. For a strong coupling between the vibrational mode and the 
thermal bath, vibrational excitation is governed by the formation of a polaronic state. 

Naturally, the molecular junction responds to an increased coupling to a "cold" thermal 
bath by decreased levels of current-induced vibrational excitation. However, as can be seen 



by inspection of Eqs. (54), vibrational excitation is not only a result of inelastic transport 
processes, but also stems from the population of the electronic states, that is the formation 
of a polaronic state [88]. Such polaron- formation leads to a finite vibrational excitation 
even in the limit of strong mode-bath coupling (i. Since the electronic level is almost fully 
populated for positive, but almost unoccupied for negative bias voltages, one observes a 
higher vibrational excitation for positive bias voltages than for negative bias voltages, if the 
mode-bath coupling strength d exceeds a value of 0.02 eV. 

We finally conclude that in an asymmetric molecular junction the level of vibrational 
excitation can be controlled by the magnitude and the polarity of the applied bias voltage. 
It is noted that a gate voltage, which allows to align the energy of electronic states, e^, with 
respect to the Fermi-level, may facilitate a control mechanism for the ratio between the 
different levels of vibrational excitation at different bias polarities. 
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B. Mode-Selective Vibrational Excitation 



In Sec. HI A we have outlined how the level of excitation of a single vibrational mode 
can be controlled by an external bias voltage $. In this section, we extend this concept 
to selective excitation of specific vibrational modes in a junction with multiple vibrational 
degrees of freedom. In particular, we show that modes with higher frequencies can be 
stronger excited than low-frequency modes, which corresponds to a "non-statistical" distri- 
bution of vibrational energy. A minimal model for two vibrational modes (model A), which 
demonstrates such mode-selective vibrational excitation, was recently introduced [31] . It in- 
volves two electronic states, where each state is coupled to one of the vibrational modes and 
asymmetrically to the leads. Thereby, the asymmetry in the coupling to the leads reflects 
an inherent asymmetry of the contacted molecule. In this section, we review and extend 
our earlier study of model A, taking into account intra-molecular correlations, in particular 
off-diagonal electronic-vibrational coupling, X V)m ^ 5 vm , and electron-electron interactions, 
U m , n 7^ 0. Moreover, we consider a different generic realization of an asymmetric molecular 
junction exhibiting MSVE, model B. Model B also comprises two vibrational modes and two 
electronic states, asymmetrically coupled to leads, but in contrast to model A, the asymme- 
try in the molecule- lead coupling is not a result of an intrinsic asymmetry of the molecule, 
but rather stems from an asymmetry in the electrodes. This corresponds for example to 
an STM setup, where the molecule bridging the gap between the two electrodes is typically 
much stronger coupled to the substrate than to the STM tip. These two scenarios, where 
MSVE can be controlled by an external bias voltage, are schematically depicted in Fig. |8j 
Respective model parameters are detailed in Table [TTJ Note that an asymmetric molecule- 
lead as well as electronic-vibrational coupling is necessary to observe MSVE in both model 
systems (cf. the discussion of Fig. |6|. 

1. The Basic Phenomenon 

First, we discuss results, where we do not account for a coupling between the vibrational 
modes and a thermal bath (d = ( 2 — 0), nor for off-diagonal electronic-vibrational coupling 
(A„, m ~ $v,m) or electron-electron interactions (U 1>2 = 0). Current- volt age characteristics 
and the corresponding population of the electronic states are shown in Figs. [9] and 11 for 
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FIG. 8: (Color online) Two generic model systems for a molecular junction exhibiting MSVE. Panel 
(a) depicts a molecular junction, where two electronic states are located above the Fermi-level of 
the junction, while panel (b) shows a molecular junction that involves an electronic state above 
and another state below the Fermi-level. Strong (weak) coupling of the electronic states to the 
leads is represented by large (small) blue dots. 



TABLE II: Model Parameters (Energy values are given in eV, K G {L,R}, a 6 {1,2}) 
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model A and B, respectively. The corresponding levels of vibrational excitation are depicted 



in Figs. 10 and 12 In both models, the excitation of the two normal modes is calculated with 
respect to the neutral molecule, in which the two electronic states (LUMO and LUMO+1) 
are unoccupied. Notice, however, that while in model A these states become occupied only 
for non-zero bias, in model B the molecule (LUMO) is charged (and therefore to some 
extent vibrationally excited) already at zero bias. Since X U)Tn oc 5 U:Tn and £7i j2 = 0, the two 
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subsystems, the one comprising state 1 and mode 1, and the other consisting of state 2 and 
mode 2, are interrelated only by the coupling of the two electronic states to the leads. This 
coupling, however, does not induce strong correlations between the two subsystems, because 
the electronic states are non-degenerate, i.e., \e 2 — e%\ > T. The transport characteristics 
of model A and B can therefore, in principle, be understood by the arguments given in 



Sec. HI A for a single electronic state and a single vibrational mode. In particular, the 



agreement of the results obtained with NEGF and ME, which we already found in Sec. 



HI A , is maintained. 

While for both models the current-voltage characteristics is almost anti-symmetric with 
respect to the applied bias voltage $, the electronic population and average levels of vi- 
brational excitation exhibit strong asymmetric behavior. In particular, for negative bias 
voltages mode 1 shows a much higher level of vibrational excitation than mode 2. For posi- 
tive bias voltages, however, the distribution of vibrational energy is reversed and mode 2 is 
higher excited than mode 1, despite the fact that Q 2 > &>i- 

These results, where no intra-molecular interactions are considered, are in line with the 
interpretation and the analysis for cooling of vibrational modes by electron-hole pair creation 



processes (cf. Sec. Ill A). In particular, it is sufficient to consider the asymmetry in the 
coupling of each electronic state to the two leads, and the energy gap between each electronic 
state and the chemical potential of the two electrodes in order to assess which of the two 
vibrations is more effectively excited. In a realistic model of a molecular junction, however, 
correlations need to be taken into account. To this end, we analyze MSVE in the next 
three sections in terms of off-diagonal electronic-vibrational coupling, X UyTn ^ 6 v , m , electron- 
electron interactions, ^ 0, and in the presence of efficient cooling by coupling to a cold 
nuclear bath, £i ^ and £2 7^ 0. 

2. MSVE in the Presence of Off-Diagonal Electronic- Vibrational Coupling 

As shown above, the MSVE phenomenon depends predominantly on the efficiency of 
cooling by electron-hole pair creation processes. This efficiency can be selectively controlled 
by the external bias voltage due to the asymmetry not only in the molecule-lead coupling, 
but also in the electronic-vibrational coupling. The latter is most pronounced when each 



mode is coupled exclusively to a different electronic state, i.e., X In Figs. \L3\ and 
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FIG. 9: (Color online) Current-voltage characteristics and the respective population of the elec- 
tronic states for the model system depicted by Fig. [8^i (model A). Solid (dashed dotted) lines are 
obtained employing the ME (NEGF) methodology. The population characteristics of state 1 and 2 
are depicted by the blue and the red line, respectively. While the current is almost anti-symmetric 
with respect to the bias voltage <&, the population of the electronic levels reflects the asymmetry 
in the coupling of the two states to the left and the right lead. 
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FIG. 10: (Color online) Average levels of vibrational excitation for the two modes in model A 
as a function of the applied bias voltage <3?. The blue and the red line depict the excitation 
characteristics of mode 1 and mode 2, respectively. Despite the different frequencies of the modes, 
mode 1 is much higher excited than mode 2 for negative bias voltages, while it is less excited 
than mode 2 for positive bias voltages. The excitation of the two modes can thus be selectively 
controlled by the external bias voltage $ (MSVE). 



14, for model A and B, respectively, the vibrational excitation of the two modes is shown for 
increasing off-diagonal electronic-vibrational coupling: X 12 = and A 2 ,i = aA 2 ,2, where 
a = describes the absence of off-diagonal vibronic coupling, while for a = 1 off-diagonal 
coupling is as strong as the diagonal one. Thereby, we use a fixed bias voltage $ = +2 V for 



Figs. |13ft and [14ft, and $ = — 2 V for Figs. [13p and [T4p. The results demonstrate that off- 
diagonal electronic-vibrational coupling tends to decrease MSVE for these model molecular 
junctions, as might have been anticipated. MSVE, however, remains significant for a broad 
range of coupling strengths a. A more detailed analysis rationalizes the trends in each case. 



In model A, at positive bias (Fig. 13 a) and for a = 0, cooling by electron-hole pair cre- 
ation is more effective via the state that is coupled more strongly to the left electrode (state 
1). Therefore, the mode coupled to this state, that is mode 1, is more effectively cooled. As 
a increases, mode 2 becomes coupled to state 1, and cooling by electron- hole pair creation 
becomes effective also for this mode. While the level of excitation of mode 1 is thus almost 
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FIG. 11: (Color online) Current-voltage characteristics and the respective population the electronic 
states for the model system depicted by Fig.[8]3 (model B). Solid (dashed dotted) lines are obtained 
employing the ME (NEGF) methodology. The population characteristics of state 1 and 2 are 
depicted by the blue and the red line, respectively. While the current is almost anti-symmetric 
with respect to the bias voltage the population of the electronic levels reflects the asymmetry 
in the coupling of the molecule to the left and the right lead. 
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FIG. 12: (Color online) Average levels of vibrational excitation for the two modes in model B as a 
function of the applied bias voltage <£. The blue and the red line depict the excitation characteristics 
of mode 1 and mode 2, respectively. Despite the different frequencies of the two modes, mode 1 
is much higher excited than mode 2 for negative bias voltages, while it is less excited than mode 
2 for positive bias voltages. The excitation of the two modes can thus be selectively controlled by 
the external bias voltage <J? (MSVE). 

the same for all values of a, the one of mode 2 decreases. Similar arguments hold for negative 
bias voltages (Fig. [Lsfc), where electron-hole pair creation via state 2 at the right electrode 
is the dominant cooling mechanism. 

Notice that off-diagonal coupling also involves a change in the nuclear reorganization energy 
of each electronic state. This is particularly pronounced for transport and pair creation 
processes involving the di-anionic states (or a doubly occupied molecular bridge, cf. Fig. 



15), which reorganization energies also involve vibrationally induced electron-electron inter- 



actions (Ui t 2 ~ 0.25 eV for a = 1, see Sec. II C). Since state 1 (2) is almost fully occupied 



for $ = 2V ($ = — 2V), processes involving state 2 (1) are dominated by the di-anionic 
resonance at 62(1) + Ui,2- Increasing a shifts this resonance to significantly lower energies. 
This shift manifests itself in the kink observed in the vibrational excitation of mode 2 (1) at 
a = 0.4 for $ = 2 V (at a = 0.2 for $ = — 2 V), indicating the suppression of electron-hole 



pair creation processes (cf. Fig. 15 1) as the respective resonance is shifted further away from 
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the chemical potential in the left (right) electrode. Note that such kinks are less pronounced 
for larger values of a, where cooling by electron-hole pair creation processes occurs for each 
mode via both electronic states such that the closure of one of these cooling channels is less 
significant. 



Similar trends are observed for model B (cf. Fig. 14). For negative bias and a = 0, mode 2 
is more effectively cooled due to the strong coupling between state 2 and the left electrode. 
Increasing a, mode 1 becomes coupled to that state as well, leading to a suppression of 
vibrational excitation also for this mode. Notice that di-anionic resonances are less important 
in this case as the two states are almost unoccupied. Therefore, kinks associated with the 
reorganization energy of the electronic levels are also less pronounced. For positive bias, 
however, the two electronic states are almost fully occupied, and therefore, transport and 
pair creation processes do occur predominantly by the di-anionic resonances. Since these 
resonances are shifted to lower energies with increasing a, state 1 is effectively located 
further away from the chemical potential in the left electrode and state 2 closer to the 
one in the right electrode. This results in less (more) efficient cooling by electron-hole 
pair creation processes, and respectively, in an increased (decreased) level of vibrational 
excitation. The latter trends bring the excitation levels of the two modes to similar values 
already for a ~ 0.2 — 0.5, which suppresses MSVE for this model. We note however, that the 
'Coulomb-like' attraction term (U), which dominates the suppression of MSVE, is typically 
compensated by repulsive electron-electron interactions, which, however, are not accounted 
for in the present model. 

3. MSVE in the Presence of Electron- Electron Interactions 

A comparison of the electronic populations, shown in Figs. [9Jd andfTTjj, and the associated 



levels of vibrational excitation, given in Figs. 10 and 12, shows that these quantities are 
strongly correlated. Electron-electron interactions, U m , n ^ 0, can strongly influence the 
electronic population of the different states [35j H22H125] and thus the degree of MSVE in 
an asymmetric molecular junction. We study the effect of such inter-state correlations on 



MSVE in Figs. [16] and 17, where the levels of excitation for the two vibrational modes in 
model A and B are plotted as functions of the electron-electron interaction strength = U, 
using fixed bias voltages $ = ±2 V. 
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FIG. 13: (Color online) Average vibrational excitation of the two vibrational modes in model A as 
a function of the off-diagonal electronic- vibrational coupling strength a. The top and the bottom 
plots correspond to a fixed bias voltage of <1> = ±2V, respectively, as illustrated in the insets. 
Although off-diagonal electronic- vibrational coupling distributes current-induced excitation among 
the vibrational modes, MSVE occurs for a broad range of coupling strengths a. 



In model A (Fig. 16), the dependence of the two vibrational mode-excitations on U is 
nearly reversed when the bias is reversed, demonstrating MSVE for most values of U. The 
particular levels of vibrational excitation and the direction of MSVE reflect the detailed 
asymmetry of the molecular junction. For example, for $ = 2 V the excitation of mode 1 
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FIG. 14: (Color online) Average vibrational excitation of the two vibrational modes in model B as 
a function of the off-diagonal electronic- vibrational coupling strength a. The top and the bottom 
plots correspond to a fixed bias voltage of = ±2V, respectively, as illustrated in the insets. 
Although off-diagonal electronic- vibrational coupling distributes current-induced excitation among 
the vibrational modes, MSVE occurs for a broad range of coupling strengths a. 

is nearly independent on U, while the one of mode 2 shows a strong non-monotonic depen- 
dence. Since mode 1 is coupled to state 1, and because state 2 is almost unoccupied in this 
regime, mode 1 is efficiently cooled by electron-hole pair creation processes with respect to 
the left electrode, regardless of the electron-electron interaction strength U. In contrast, 
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FIG. 15: (Color online) Schematic representation of Coulomb-assisted transport processes (Panel 
(a)) and Coulomb-assisted electron- hole pair creation processes (Panels (b) and (c)). The dashed 
line represents the position of the di-anionic (doubly occupied) resonance, where state 2 is consid- 
ered occupied. Similar processes are available at energies e 2 + U. The position of the two electronic 
states with respect to the chemical potentials in the left and the right lead, e l7 / 2 — //l/R; an d the 
positions of the di-anionic resonances, e ly / 2 + U — /U L / R , determine the efficiency of heating and 
cooling processes. 

mode 2 is coupled to state 2 that is only weakly coupled to the left lead, and therefore, 
not efficiently cooled by electron-hole pair creation processes for U — 0. However, as U 
increases, electron tunneling with respect to state 2 involves an increasingly larger charging 
energy, I2+U, since state 1 is almost fully occupied. This brings the electronic energy in this 
transport channel first closer to and then further away from the chemical potential in the 
left electrode. Cooling of mode 2 by electron-hole pair creation processes (Coulomb-assisted 



electron-hole pair creation processes as depicted in Figs. 15 d and 15 2) is thus first enhanced 
and then suppressed as U increases, resulting in the observed non-monotonic level of vibra- 
tional excitation. Notice that the pronounced cooling due to Coulomb-assisted electron-hole 
pair creation processes, in the intermediate regime of electron-electron interaction strengths, 



U, reverses the direction of MSVE with respect to the U = case (cf. Fig. 16) 



In model B (Fig. 17) the asymmetry in the coupling of the molecular bridge to the 
electrodes leads to very different dependencies of the vibrational mode-excitations on U for 
different polarities of the applied bias voltage. The resulting dependence of MSVE on U is 
non-trivial, ranging from enhancement to suppression of the effect with respect to U = 0. 
At negative bias voltages the two electronic states remain nearly unoccupied, so that the 
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FIG. 16: (Color online) Vibrational excitation as a function of the electronic interaction strength 
U in model A. Blue (red) lines refer to the average excitation number of mode 1 (2). The top and 
the bottom plots correspond to a fixed bias voltage of = ±2V, respectively. The efficiency of 
cooling mode 2 (1) by electron- hole pair creation processes with respect to the di- anionic (doubly 
occupied) state varies with the electron-electron interaction strengths U for positive (negative) bias 
voltages. This leads to a inversion of MSVE for intermediate values of U, where the energy of the 
corresponding di-anionic state, e2 + U (ei + U), is close to the chemical potential in the left (right) 
lead. 
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FIG. 17: (Color online) Vibrational excitation as a function of the electronic interaction strength 
U in model B. Blue (red) lines refer to the average excitation number of mode 1 (2). The top and 
the bottom plots correspond to a fixed bias voltage of $ = ±2V, respectively. For negative bias 
voltages, due to the weak coupling to the right lead, both states are more or less unpopulated, 
which translates to a weak dependence of vibrational excitation on U. For positive bias voltages, 
however, double occupancy of the molecular bridge becomes important, and thus, the specific levels 
of vibrational excitation show a strong dependence on U. While MSVE thus exhibits a strong 
enhancement for weak and intermediate electron-electron interactions strengths, it is attenuated 
for higher values of U. 
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effect of electron-electron interactions is negligible. The lower excitation level of mode 2 at 
this polarity of the bias voltage reflects more efficient cooling by electron-hole pair creation 



via state 2 at the left electrode (Fig. 17), which is maintained for different values of U. 
For positive bias voltages, the two electronic states are almost fully occupied. Cooling of 
mode 1 and 2 is thus dominated by electron-hole pair creation at the left and the right leads, 
respectively, according to the proximity of the corresponding energy levels, e\ + U and I2 + U , 
to the respective chemical potentials. As U increases, electron-hole pair creation at the left 
electrode is enhanced, while pair creation with respect to the right lead is suppressed. This 
leads first to an enhancement of MSVE with respect to U — 0. At U ~ 0.4 eV, however, the 
electronic energy ei + U crosses the chemical potential in the left lead, and state 1 becomes 
discharged. At this point, cooling of mode 1 is at its maximal efficiency. Simultaneously, 
since state 1 is no longer populated ((ajai) drops from « 0.9 to w 0.1), processes involving 
state 2 occur predominantly via the anionic channel at £2, such that cooling of mode 2 by 
electron-hole pair creation processes with respect to the right electrode becomes as efficient 
as for U = 0. Accordingly, the level of excitation for mode 2 returns to its original value. 
Increasing U even further cooling of mode 1 (via the di-anionic channel, e% + U) becomes 
less efficient, and the level of excitation for mode 1 increases again. As €2 + U is closer 
to the chemical potential in the left lead, the cooling efficiency of mode 2 by electron-hole 
pair creation processes with respect to the left lead increases. The overall effect leads to 
a suppression of MSVE for U > leV in the given range of electron-electron interactions 
strengths U. Note, however, that for yet larger values of the electron-electron interaction 
strength, U > 2.5 eV, MSVE is regained and becomes approximately as pronounced as for 
U = 0. In this regime, transport and pair creation processes are dominated by the anionic 
resonances (at 61/2) so that the asymmetry in the cooling efficiency of the two modes (at 
least for the present model) is the same as for U = 0. 

4- MSVE for Strong Vibrational Relaxation 



In the discussion of Fig. [7] in Sec. HI A , we have already seen that the level of excitation of 



a single vibrational mode consists of two contributions: current-induced local heating due to 
inelastic electron transport processes (as shown in Figs. (2^,-d), and polaron-formation |88j, 
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which can be quantified by the difference 

( c I c v)h - ( c Uu)h = ^{K, m K,m' /^l){a ] m a m a} m ,a m >) R . (58) 

mm 1 

Strong vibrational relaxation results in a strong suppression of current-induced vibrational 
excitation, especially if the time-scale for vibrational relaxation is much shorter than the 
time-scale between two consecutive transport events. The contribution due to polaron forma- 
tion, however, stems from the steady-state population of the electronic levels in a molecular 
junction. Hence, for model A, MSVE occurs in the presence of strong vibrational relaxation 
[31], since the population of the electronic levels can be selectively controlled by the external 
bias voltage $ (cf. Fig. ^p)- In model B, however, both states are either fully populated 
or empty such that strong vibrational relaxation is likely to hinder MSVE for this model 
system. 

IV. CONCLUSION 

In this work we have studied and analyzed transport characteristics of single-molecule 
junctions, focusing on the excitation of specific molecular vibrational modes. In particular, 
we have shown that the level of excitation of specific modes can be controlled by the polarity 
and the magnitude of an external bias voltage. Thereby, high-frequency modes (typically 
associated with strong chemical bonds) can be higher excited than low-frequency modes, 
which translates to a "non-statistical" distribution of energy among the vibrational modes. 
We refer to this phenomenon as mode-selective vibrational excitation. 

Our main findings are summarized below: 

1) The importance of cooling by electron-hole pair creation 

Our analysis shows that cooling of the vibrational modes in a molecular junction by 
electron-hole pair creation processes is crucial to understand the extent of the MSVE 
phenomenon. In particular, since the efficiency of these processes is sensitive to the 
position of the chemical potentials in the leads, the levels of vibrational excitation in 
the molecule can be controlled by an external bias voltage. Considering a molecule 
with multiple vibrational modes and typical asymmetries in the vibronic as well as 
molecule-lead couplings, the level of excitation of specific vibrational modes can thus 
be tuned by the polarity and the magnitude of the external bias voltage. 
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2) The role of asymmetry and intra-molecular interactions 

Our studies suggest that MSVE is a rather general phenomenon and is likely to be 
observed experimentally. The required asymmetry in the electronic interaction be- 
tween different molecular states and the leads may be due to an inherent asymmetric 
molecular structure (model A) or stem from an inherent difference between the two 
electrodes, as e.g. in STM experiments (model B). Intra-molecular interactions, e.g. 
due to off diagonal electronic- vibrational coupling or electron-electron repulsion tend to 
redistribute the excitation energy between the different modes, and thus work against 
MSVE. However, having analyzed a broad range of parameters, we found the MSVE 
phenomenon to prevail even in the presence of such interactions. 

3) The importance of off-resonant processes: Comparing ME to NEGF calculations. 
Our numerical studies of generic models of molecular junctions were based on two 
complementary theoretical methods: a nonequilibrium Green's function approach [32T 
134"] 138] and a master equation approach [561 EH EE] . Both approaches are based on 
a second-order expansion in the coupling of the molecular bridge to the leads. While 
the NEGF method also accounts for higher-order effects, the ME approach describes 
only resonant electron tunneling processes. However, intra-molecular interactions, 
either due to electronic-vibrational coupling or electron-electron interactions, can be 
described by the ME method without invoking further approximations, while these 
interactions are described by our NEGF approach approximately in terms of a non- 
perturbative scheme. Although the results obtained by both methodologies agree in 
most cases reasonably well, further insights into the relevant mechanisms can be gained 
when the results exhibit differences. Thus, for example, the importance of off-resonant 
electron-hole pair creation processes for local cooling [331 E3 I126H128| of vibrational 
modes in the high-bias regime could be revealed. 

We end by noting that this work considered only generic models to study the basic mech- 
anisms and prerequisites of bias-controlled MSVE. The identification of specific molecules 
that exhibit MSVE requires transport studies based on first-principles electronic structure 
calculations [80 ] I110H1121 1112H115] . This will be the subject of future work. Another inter- 
esting extension concerns the external control mechanism for MSVE. In the present work, 
we have considered an external bias voltage as the means to control MSVE. A gate elec- 
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trode [20j H16H118] may provide another tool for addressing a molecular junction with an 
electric field, and thus, may also be used to control vibrational excitation. Finally, in the 
context of mode-selective chemistry, studies of MSVE in single-molecule junctions may pave 
the way to control chemical processes in molecules adsorbed on surfaces. For example, an 
additional electrode that provides an external potential bias may induce catalytic reactions 
in a selective manner. 
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